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We concisely review the recent evolution in the study of parafermions—exotic emergent excitations 
that generalize Majorana fermions and similarly underpin a host of novel phenomena. First we illus¬ 
trate the intimate connection between Z 3 -symmetric “spin” chains and one-dimensional parafermion 
lattice models, highlighting how the latter host a topological phase featuring protected edge zero 
modes. We then tour several blueprints for the laboratory realization of parafermion zero modes— 
focusing on quantum Hall/superconductor hybrids, quantum Hall bilayers, and two-dimensional 
topological insulators—and describe striking experimental fingerprints that they provide. Finally, 
we discuss how coupled parafermion arrays in quantum Hall architectures yield topological phases 
that potentially furnish hardware for a universal, intrinsically decoherence-free quantum computer. 


I. INTRODUCTION 

The theoretical and experimental search for topologi¬ 
cally nontrivial condensed-matter systems with fraction¬ 
alized excitations continues to attract widespread at¬ 
tention. One reason for the excitement is fundamen¬ 
tal, as fractionalization brilliantly illustrates the rich 
emergent long-distance behavior that many-body sys¬ 
tems can exhibit. Another reason is more practical: 
certain ‘non-Abelian’ fractionalized excitations form the 
bedrock of topological quantum computers that tantaliz- 
ingly promise inherent immunity against errors 1 2 . 

Although fractionalized topological phases appear 
quite exotic, several fundamental examples of such states 
share a profound relationship with the very well known 
Ising quantum spin chain. One facet of this relationship 
arises from the Jordan-Wigner transformation, a non¬ 
local map between the bosonic Ising spins and fermions. 
While this mapping has long been known to provide a 
useful tool for solving the two-dimensional classical Ising 
model and its quantum-chain limit 3 -®, great insights arise 
by treating the fermions themselves as physical degrees 
of freedom. With the latter viewpoint the quantum 
Ising chain corresponds to an electronic model for a one¬ 
dimensional (ID) p -wave superconductor. 

Remarkably, conventional spin ordering for the Ising 
chain yields a topological superconducting phase in 
fermionic languages The degenerate ‘up’ and ‘down’ or¬ 
dered spin configurations of the former translate into 
even- and odd-fermion-parity superconducting states. 
However, the degeneracy in the latter is resilient; a local 
magnetic field that readily splits the energy of the spin 
polarizations becomes non-local in fermionic language. 
The p-wave superconductor thus exhibits a ground-state 
degeneracy immune to physical perturbations—one of 
the hallmarks of a topological phase. 

Because of their inherent robustness, the supercon¬ 
ducting ground states can encode fault-tolerant ‘topo¬ 
logical qubits’. The qubits arise from Majorana zero 
modes that bind to the superconductor’s edges.® (Zero 


modes in the Ising spin chain with open boundary con¬ 
ditions were found long agcP, but their topological sig¬ 
nificance not appreciated until Kitaev’s seminal work.) 
More interestingly still, the superconductor supports the 
remarkable phenomenon of non-Abelian exchange statis¬ 
tics; that is, swapping pairs of zero modes rotates the sys¬ 
tem’s quantum state—and hence the topological qubits— 
in exquisitely precise waysP® These highly attractive 
properties have inspired enormous activity; recent re¬ 
views of this fast-moving field appear in Refs. [T01IT51 


The close connection between spin systems and topo¬ 
logical phases extends in interesting and often sur¬ 
prising ways to ‘clock models’ with Z n symmetr^l 
These models admit a non-local mapping to chains of 
parafermions that comprise Z n generalizations of Majo¬ 
rana operators 17 . Our goal here is to review the fan¬ 
tastically rich physics that has emerged through recent 
studies of these more exotic objects. In many ways the 
story parallels the Ising/Majorana problem highlighted 
above, but certainly not trivially so; throughout we will 
encounter numerous features unique to parafermions. 


An immediately glaring distinction is that all of 
condensed matter arises from conventional bosons and 
fermions, so that explorations of parafermions might 
naively seem purely academic. We stress that this is 
far from the case. A substantial fraction of our discus¬ 
sion will survey plausible experimental routes to trap¬ 
ping parafermionic excitations even in well understood, 
presently available systems. Although more challenging 
to implement, these platforms give rise to topological 
qubits that are both better protected against environ¬ 
mental noise and allow for richer fault-tolerant qubit ro¬ 
tations compared to Majorana-based architectures. In 
fact we will eventually describe methods for leveraging 
systems with parafermions to synthesize fully universal 
topological quantum computing hardware! 
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II. PARAFERMIONS IN CHAINS 

In this section we review the basics of spin chains 
with Z 3 symmetry and their description in terms of 
parafermions. While some of the analysis straightfor¬ 
wardly generalizes the Z 2 Ising/Majorana case, we high¬ 
light several facets not present there. Our discussion both 
lays the groundwork for the novel experimental propos¬ 
als discussed in subsequent sections and also illuminates 
several problems of more formal interest. 


A. Lattice parafermions 

Clock models provide a natural Z n -symmetric gener¬ 
alization of the Ising model. They exhibit similar or¬ 
dered and disordered phases, and also admit a non-local 
representation wherein symmetry-breaking states map to 
topological phases supporting local zero modes 16 . Addi¬ 
tionally, for n > 2 a variety of interesting behavior not 
present in the n = 2 Ising case emerges, especially when 
the interactions are chiral 18 21 . Th e co rresponding topo¬ 
logical phases have been classified 22 1 23 1 . 

We focus primarily on the ID three-state clock model 
since it contains almost all the physics of interest. The 
Hilbert space consists of a three-state “spin” on each site. 
The usual Pauli matrices a z and a x respectively gener¬ 
alize to a and r, defined by 



Then the “clock” operator cr a = l 0 --- 0 l 0 cr 0 l... 
measures the spin on site a, while the “shift” operator 
T a = l0---0l0r0l... winds its value by uj. These 
operators obey the algebra 

= 1 , li = 1 , cr a T a = UT a o a ( 2 ) 

and commute on different sites, e.g., o a T\) — Tbcr a for 
a^b. 

The most general nearest-neighbor Hamiltonian for the 
three-state clock model on an open L-site chain is 

h = - j J2 + b.c.) PM + fj.c.). 

a=1 a=1 

(3) 

Throughout we assume purely real couplings J, /, p and 
0, with / and J non-negative. The Hamiltonian preserves 
several symmetries; see Ref. [S3] for a detailed discussion. 
As a particularly important example, there is a global Z 3 
symmetry corresponding to a uniform shift of all spins 
generated by Q = Yl a r which satisfies Q 3 = 1 . The 
Hamiltonian also possesses duality symmetries when / = 
J and (p = 6 that strongly constrain the phase diagram. 


The ferromagnetic case p = 0mod27r/3—usually 
called the three-state Potts model—has been exception¬ 
ally well-studied from many viewpoints^ some recent ef¬ 
forts motivated by the connection to topological physics 
include Refs. [22j [26j and [27j The model exhibits a self¬ 
dual critical phase transition at / = J separating ordered 
(/ < J) and disordered (/ > J) phases. In the ordered 
phase the system admits three degenerate ground states 
characterized by (cr a ) ^ 0. Acting with the generator 
Q cycles through these ground states, indicating sponta¬ 
neous Z 3 symmetry breaking. The disordered phase, by 
contrast, yields a unique, symmetric ground state with 
(Co) = 0. 

When <j> 7 ^ 0 mod 7 r/ 3 , the clock model model is chi- 
ral, breaking spatial parity and time-reversal symme¬ 
tries. The chiral case exhibits a much more intricate 
phase structure than the ferromagnetic limit. Some in¬ 
formation can be obtained by exploiting integrability on 
the two-parameter manifold J cos(30) = /cos(30), where 
for example an incommensurate phase with ground-state 
level crossings appear^. A systematic numerical anal¬ 
ysis was (finally!) performed recently, finding accurate 
phase boundaries between ordered, trivially disordered, 
and incommensurate phased^. Interestingly, there is evi¬ 
dence that incommensurability occurs along the self-dual 
/ = J, <p = 0 line all the way up to the ferromagnetic 
limit. 

Broken symmetry states of the clock model deeply re¬ 
late to exotic topologically nontrivial phases of matter—a 
remarkable correspondence that becomes manifest upon 
exploiting the Fradkin-Kadanoff mapping 17 that rewrites 
the Hamiltonian in terms of parafermions. As with Ma- 
jorana fermions in the Ising chain, parafermions are de¬ 
fined via products of ‘order operators’ cr a and ‘disorder 
operators’ built from r-strings that create kinks by wind¬ 
ing a macroscopic set of spins. Precisely, one set of 
parafermion operators reads 



0^2a-l 

11 

V 

aV 

+ 

+ 

to 



a 

= (°'i) T a t +l T a t +2 • • • . 

(4) 

which obey 




a b 

= 1 , 

a c = cd sgn(c - 6) a c a & . 

(5) 


Notice that the parafermions do not commute at different 
sites; this important property reflects their non-locality 
in clock variables. Equally importantly, they do not an¬ 
ticommute either, unlike fermions. Instead exchanging 
parafermion operators yields a complex phase of uo —a 
situation reminiscent of anyons. 

In terms of parafermions, the Hamiltonian is 

H = - Je* £ u 2 al + 1 a b - fe ie £ M b+1 a b + f».c. 

b even b odd 

(6) 

Even though H is bilinear in parafermions, the factor 
of in the algebra ^ means that the model cannot 
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be solved simply by Fourier transformation (contrary to 
the Ising case). It is worth noting however that the en¬ 
tire spectrum of the non-Hermitian Hamiltonian found 
by omitting ‘-hlpc.’ in Eq. © is in a precise sense that of 
“free” parafermion d 29 * 30 -! 

At this point an interesting question arises: exactly 
how is the physics of the ordered and disordered clock- 
model phases encoded in the parafermion representa¬ 
tion? The answer directly relates to the existence of 
parafermion zero modes, which play an essential role 
throughout this review. 


B. Edge zero modes 

One important signal of a topological phase in ID is 
the existence of robust edge zero modes, the analogue 
of gapless modes in higher dimensions. To be precise, 
consider a system with a discrete symmetry generated 
by an operator Q. We will say that the system pos¬ 
sesses a ‘strong zero mode’ if there exists an operator that 
commutes with the microscopic Hamiltonian up to terms 
falling off exponentially with system size, but does not 
commute with Q. The energy levels for each eigenvalue 
of Q must then be identical (modulo exponentially small 
corrections), resulting in degeneracies across the entire 
spectrum. Strong edge zero modes have the same prop¬ 
erty but additionally localize with exponential fall-off to 
the boundaries of an open system. As a famous example, 
the topological phase in the Major ana-chain representa¬ 
tion of the Ising model supports strong edge Major ana 
zero modes, which one can readily construct explicitly 
and exactly!^ 

We emphasize that strong edge zero modes impose an 
extraordinarily strong constraint on the spectrum that 
propagates to arbitrarily high energies. In glaring con¬ 
trast, a ID topological phase in general requires only 
ground-state degeneracy—a much weaker condition. It 
is thus useful to introduce the notion of a ‘weak edge 
zero mode’ that only commutes with the Hamiltonian 
projected into a low-energy subspace, thus encoding the 
minimal degeneracy necessary for a topological phase. 
While harder to find, strong zero modes do offer poten¬ 
tial advantages, since one in principle could exploit the 
excited-state degeneracies to build a fault-tolerant qubit 
that works even at high temperatures.^ 

With these preliminaries in place, let us return to the 
parafermion representation of the clock chain. To estab¬ 
lish a baseline we first examine two limits, illustrated 
in Figs, [lja) and (b), where the model is trivially solv¬ 
able. In the figure circles within each box represent the 
two parafermions o^a-i and ot 2 a associated with a single 
clock site a; recall Eqs. ©• 

In the maximally disordered limit J = 0, spins on dif¬ 
ferent sites simply decouple, and so parafermions interact 
only with their partner on the same clock site as indicated 
by the bonds in Fig. [lja). The parafermion chain then 
generically realizes a gapped trivial phase with a unique 


(a) J = 0, / > 0 (trivial) 



(b) J > 0, / = 0 (topological) 



(c) J = f > 0, (j) = 6 = 0 (critical) 

- > - 



FIG. 1. Schematic representation of the parafermion chain 
[Eq. deep in the (a) trivial and (b) topological phases, 
and (cj at the intervening critical point. 


ground state, consistent with the physics in the spin rep¬ 
resentation. 

When / = 0 the model maximally orders in spin lan¬ 
guage; accordingly, parafermions on adjacent clock sites 
strongly couple as represented by the bonds from one 
box to the next in Fig. 0 b). Here too the bulk of 
the parafermion chain generically acquires a gap, though 
now the edges behave more interestingly. In particu¬ 
lar, a novel feature at / = 0 is that the parafermions 
aq and a^z, on the ends disappear from the Hamilto¬ 
nian and thus commute with it. They do not, however, 
commute with the Z 3 symmetry generator Q. The sys¬ 
tem therefore exhibits a pair of strong edge parafermion 
zero modes ol\,ol 2 l that signify the onset of a topologi¬ 
cal phase. When acting on any mutual eigenstate of H 
and Q, the zero modes yield a distinct eigenstate with the 
same energy. (The states must differ because if q denotes 
the initial eigenvalue of Q, acting with aq or a^z, cycles 
the eigenvalue to uoq.) Thus the ground states and all ex¬ 
cited states assume at least threefold degeneracy in the 
/ = 0 limit. While the degeneracy has a topological ori¬ 
gin in the parafermion representation, in spin language it 
arises from the three sectors of symmetry-breaking states 
with (a a ) 7 ^ 0 and their excitations. 

The fate of the strong parafermion zero modes away 
from this trivial limit raises an important question: are 
there modified operators oq + • • • and OL 2 L + • • • that 
commute with the Hamiltonian even when / 7 ^ 0? A 
natural guess is that strong zero modes persist pro¬ 
vided the bulk gap remains finite so that the topologi¬ 
cal phase is intact; this intuition indeed proves correct in 
the Ising/Majorana case. The interesting answer here, 
however, is that the strong zero modes require chirality 
(i.e., (j) 7 ^ 0 mod 7 r/ 3 ) to survive® This result is particu¬ 
larly surprising given that the topological phase is most 
robust (i.e., has the largest gap) in the ferromagnetic 
limit and, more broadly, always remains for sufficiently 



























4 


small / as long as the system orders ferromagnetically 
in spin language . 22 27 32 . (Any topological significance of 
the incommensurate phase is a so-far unexplored issue.) 
Instead, the topological phase in the non-chiral ferromag¬ 
netic 0 = 0 limit generically supports only weak edge 
parafermion zero modes of the form aq ? 2 l+* • • that nicely 
shift between three degenerate ground states but are not 
especially meaningful when acting on excited states P3 

There are several ways of seeing why strong 
parafermion zero modes require chirality. One method 
utilizes an iterative procedural Namely, split the 
Hamiltonian in Eq. (J6| into two pieces H = JV + fF 
in the obvious way, so that [oq, V] =0 (and similarly for 
ql2l)- Because [iif, oq] = f [F, oq] is non-vanishing, the 
first step is to find some X satisfying [F, aq] = [V, X}. 
If such an X exists, then [H,a i — (f /J) X] is of order 
/ 2 . One can attempt to repeat this procedure to con¬ 
struct strong zero modes as a power series in f /J. This 
exercise is straightforward for the Ising/Major ana cas^ 
however, one finds 16 1 for the Z 3 parafermion chain that 
X is proportional to 1/sin(30). Thus the iteration im¬ 
mediately breaks down for the non-chiral ferromagnetic 
limit—despite the robust topological bulk gap—implying 
the demise of strong edge zero modes at arbitrarily small 
but finite f/J. Conversely, Ref. HBl proved that the it¬ 
erative procedure can be implemented at all orders for 
sin(30) 7 ^ 0. Strong edge zero modes therefore exist for 
sufficiently small //[Jsin(30)], although an explicit all¬ 
orders expression is not known. As in the Ising/Majorana 
case, the uniform couplings considered here are not ac¬ 
tually necessary; the above arguments still work when / 
and J vary spatially. 

Perturbation theory in / provides a complementary 
and illuminating way of understanding this curious 
resultP 2 ! First, we verify that the ground states in the 
topological phase remain exponentially degenerate for 
sufficiently small //J, as required for weak parafermion 
zero modes and hence topological behavior to exist. For 
the following it proves convenient to revert back to the 
equivalent clock-model representation and label the three 
spin orientations at each site by A, B , and C. At / = 0 
and |0| < 7 r /3 the system exhibits three ground states 
given simply by |AAA ...), |BBB ...), and |CCC ...). 
Since the / perturbation winds single spins, its action 
on, say, the A ground state gives a sum over excited 
states of the form fe l9 (\BAA ...) + | ABAA ...) + ...) + 
fe~ l0 (\CAA ...) + \ACAA all having an ex¬ 
citation energy of order J. The first perturbative cor¬ 
rection to the ground-state energies then arises at order 
/ 2 /J, but crucially, Z 3 symmetry requires it be identical 
for all three ground states. The energy splitting between 
different ground states occurs only at order L in per¬ 
turbation theory because mixing ground states requires 
shifting all L spins. This yields a splitting proportional to 
(//J) L , indeed exponentially small in system size. Such 
degeneracy is all that is necessary for the persistence of 
a topological phase in the parafermion chain. 

Qualitatively different behavior appears for the ex¬ 


citations. Consider the ferromagnetic 0 = 0 limit. 
The lowest-lying excitations arise from six flavors 
of ‘one-kink’ states of the form | ... AAABBB ...), 

| ... BBB AAA ...), etc., all exactly degenerate at / = 0. 
Winding of single spins induced by the / perturbation 
not only allows the kinks to hop, but in an open chain 
also enables efficient conversion of one flavor to another. 
More precisely, sequential actions of / terms can produce 
the following sequence: (i) An | ... AAABBB ...) kink 
moves to the right end of the system where it becomes 
| ...AAAB), (ii) winding the rightmost spin changes 
kink flavor to \ ...AAAC), and (Hi) hopping the kink 
leftward similarly allows | ACCC | BCCC ...) 

| ... BBBCCC ...). The key is that we connected states 
related by a global Z 3 transformation—precisely those 
linked by strong parafermion zero modes at / = 0 — 
without leaving the degenerate subspace. These exci¬ 
tations thus acquire power-law splitting at any finite /, 
implying the immediate demotion of strong to weak zero 
modesP^ 

In the chiral case, broken spatial parity causes the en¬ 
ergy of the intermediate | ... AAACCC ...) kink states 
to differ from that of the initial | ... AAABBB ...) state. 
An enormous energy barrier consequently develops that 
sharply suppresses the above process and drives expo¬ 
nential instead of power-law splitting for the low-lying 
excitations. Thus chirality naturally favors the revival 
of strong zero modes, consistent with the iterative con¬ 
struction above. 

Intriguingly, a sharp transition (tuned, e.g., by the chi¬ 
ral phase 0 ) appears to separate couplings that produce 
strong vs. weak edge zero modes. Detailed numerics, for 
instance, do indeed see abrupt changes between expo¬ 
nential and power-law splitting among excited states.^ 
Even more interestingly, the numerics suggest a contin¬ 
uous transition, in the sense that the splittings conform 
well to a scaling ansatz with ‘critical’ exponents. A par¬ 
ticularly striking feature worth emphasizing is that this 
transition is purely a property of excited states; indeed 
the ground states in the topological phase always exhibit 
exponential splitting and show no signs of ‘criticality’. 
In this sense the physics is reminiscent of transitions dis¬ 
cussed in the context of many-body localization 31 33 . 


C. Parafermions in conformal field theory 

Conformal field theory enters the physics of topolog¬ 
ical phases in multiple ways. It provides a way of clas¬ 
sifying all the operators and characterizing their opera¬ 
tor products—i.e., their “fusion”—in many ID quantum 
critical points. Such theories, for instance, describe con¬ 
tinuous transitions out of topological phases, including 
in the parafermion chain. They also describe the gapless 
edges of many two-dimensional topological phases. Re¬ 
markably, the same underlying mathematical structure 
remains important even away from quantum criticality. 
In fact it also characterizes the fusion and braiding of 
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anyons in the gapped 2D topological phase that we will 
construct in Sec.[lV]by coupling critical chains together. 

Parafermionic conformal field theory, introduced and 
described in depth by Zamolodchikov and Fateev 34 !, gov¬ 
erns the scaling limit of the three-state Potts ferromag¬ 
netic critical point [f = J and 9 = (j) = 0 in Eqs. © 
and ([6|] that separates the topological and trivial phases 
of the parafermion chain; see Fig. |TJc) for an illustra¬ 
tion. Not surprisingly, there are parafermionic fields— 
conventionally labeled V+ and ijjR —related to the lattice 
parafermion operators. The former left-moving field de¬ 
pends on the space-time coordinates x and t solely via the 
combination x + vt for some velocity v, while the latter 
right-mover depends solely on x — vt. Each exhibits scal¬ 
ing dimension 2/3. Like their lattice counterparts, they 
obey the operator products ifp X ~ and V’E ^ 1, 
and likewise for i/jr. 

More surprisingly, taking the continuum limit of the 
lattice operator does not give the parafermion fields 
V+ and/or 'i/jr as the le ading pieces; other more relevant 
operators mix in as welP^ 35 . It is both formally inter¬ 
esting, and essential for the construction of 2D topolog¬ 
ically ordered phases reviewed in Sec. |IV[ to find lattice 
operators that admit exclusively ^ or ^r as their domi¬ 
nant contribution in the long-distance limit, i.e., to filter 
out unwanted terms. Fortunately this is indeed possi¬ 
ble. A detailed analysis utilizing the problem’s discrete 
symmetries shows that the continuum limit of yields 
the combination 


~ ciip L + c 2 (-l) b $ , (7) 

where c\ and C 2 denote unimportant numerical constants 
and <F is the unwanted field (which carries dimension 
1/15 + 2/5 = 7/15 < 2/3). Because of the alternat¬ 
ing sign in the last term, a simple linear combination 
successfully isolates the parafermion field: 

«2a-l + «2a ~ V+ • (8) 

Finding a lattice analog of 'i/jr requires introduc¬ 
ing another—though not independent—set of lattice 
parafermion operators, 

fi2a—l = •••+2 + i(<+, 

@2 a = ■■■ + 2+1 ++ 1 ) • ( 9 ) 

Compared to Eqs. Q these operators simply have r 
strings running in the opposite direction. The cds and 
/3’s are hence related by spatial parity, which immedi¬ 
ately implies that 


scaling dimension 4/5. This Z 3 -invariant field describes 
the scaling limit of the perturbation that moves the criti¬ 
cal parafermion chain into either the topological or trivial 
phases. In notation introduced in Sec. IIB the pertur¬ 
bation corresponds to the lattice operator F — V with 
0 = (j) = 0. Fusion of the energy field takes the form 
E x E = 1 + E, which underlies the appearance of ‘Fi¬ 


bonacci anyons’ that we will encounter in Sec. IV 


III. EXPERIMENTAL BLUEPRINTS FOR 
PARAFERMION ZERO MODES 

We turn now to the crucial question of how one can 
experimentally realize parafermion zero modes. In what 
follows we focus exclusively on weak zero modes, though 
for brevity we hereafter drop the ‘weak’ qualifier. Even 
for this case the problem poses a significant challenge— 
certainly more so compared to the pursuit of their Ma- 
jorana cousins. For one, the parafermion chain Hamilto¬ 
nian constructed in Sec. Ill Al involves neither bosons nor 
fermions as elementary constituents; thus at a minimum 
parafermion zero modes require a strongly correlated host 
featuring nontrivial emergent degrees of freedom. The 
Majorana chain, by contrast, caricatures a p-wave super¬ 
conductor built from electrons and admits a simple free- 
particle description. Further subtlety awaits even after 
we abandon the comforts of free particles: rigorous classi¬ 
fications of ID bosonic and fermionic systems 3 ^^ (sub¬ 
ject to a restricted set of symmetries) capture phases sup¬ 
porting Majorana zero modes but nothing more exotic— 
even with strong interactions! 

Fortunately loopholes exist. A recurring theme of this 
section will be the use of ‘anomalous’ edge states in 2D 
systems as a cradle for ID phases that host parafermion 
zero modes. By definition such edge states can not 
appear in strictly ID media and hence fall outside of 
the aforementioned classifications. Below we survey sev¬ 
eral candidate platforms, focusing on those built upon 
presently available physical systems. We stress that other 
interesting, often very closely related proposals exist in¬ 
cluding an important early work that invokes lattice de¬ 
fects in fractional Chern insulator^! generalized rotor 
models®! and fractionalized topological insulators and 
superconductors 4 ®®! (among others). 


A. Quantum Hall/superconductor hybrids 


1. Setup and physical origin of zero modes 


ha-1 + 02a ~ i>R ( 10 ) 

as desired. 

A short overview of the parafermion conformal field 
theory’s fusion algebra, geared toward the present con¬ 
text, appears in Ref. [35j One feature worth noting is the 
Fibonacci structure inherent in the “energy” field E of 


The ‘anyonic’ commutation relation exhibited by the 
zero modes of interest [recall Eq. ©] spotlights fraction¬ 
alized 2D systems with anyons as natural parafermion 
hosts. Fractional quantum Hall phases provide many ex¬ 
perimentally relevant examples, and luckily one can trap 
parafermions even in the simplest varieties. For peda¬ 
gogical purposes we will first examine in some detail a 
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(c) 


Superconductor 



FIG. 2. Select blueprints for parafermion zero modes a± : 2 - In (a) only the charged edge mode is shown for simplicity. 


platform that closely parallels various well known Ma- 
jorana proposals 45 47 : a spin-unpolarized quantum Hall 
state at filling v = 2/3 (i.e., the 112 state) interfaced with 
an s-wave superconductor. Earlier discussions of other 
fractions—which closely inform our treatment below— 
appear in Refs. (48| and 491 

Figure [ 2 |a) sketches the setup. A trench carved into 
the v = 2/3 fluid creates a fractionalized ‘wire’ consisting 
of edge states that, importantly, support gapless avatars 
of the anyons living in the adjacent bulk. Since an identi¬ 
cal set of counterpropagating modes appear opposite the 
trench, physical perturbations can always gap them out. 
Electron tunneling across the divide provides one mecha¬ 
nism that simply fuses the two fluids together. Filling the 
trench with an ordinary s -wave superconductor allows 
a second, intuitively quite different possibility: generat¬ 
ing a gap through proximity-induced spin-singlet Cooper 
pairing of electrons from the two edges. These tunneling 
and pairing mechanisms gap the trench’s spin sector in 
the same way but yield topologically distinct ID gapped 
phases for the charge sector. Domain walls separating 
the two phases as in Fig. |2|a) bind zero-energy modes 
that, owing to the system’s fractionalization, correspond 
to Z 3 parafermions. (For a detailed analysis see Refs. 1351 
andl5UP 

This setup bears the remarkable feature that the spin- 
singlet Cooper pairs induced by the parent supercon¬ 
ductor arise more fundamentally from conglomerates of 
2e/3 charges present in the quantum Hall edges. The 
proximity effect therefore not only condenses charge- 2 e 
bosons at the trench, but crucially also catalyzes ‘anyon 
condensation®]. To make this precise, we define V^e /3 to 
be the spin-conserving quasiparticle operator that adds 
charge 2e/3 symmetrically to the trench, with e/3 on 
each end. Then in a pairing-gapped region 

2e/z) ~ e*^“ P2e/3 , (11) 


where p 2 e /3 takes on integer values. Thus bi’le/i)' which 
essentially yields the expectation value of an ordinary 
charge- 2 e Cooper pair, becomes ^ 1 . A further remark¬ 
able implication follows: although the parent supercon¬ 
ductor can only absorb paired electrons, the trench can 
subsume fractional 2e/3 charges without energy penalty! 

An analogous situation emerges in tunneling-gapped 
regions. Electron hopping—which also involves multiple 
fractional charges—catalyzes a different flavor of anyon 
condensation with 

(^dipole) ~ e i ^ Pdipole , (12) 

where V’dipoie creates a charge difference of 2e/3 across 
the trench and Pdipoie €= Tunneling-gapped regions ac¬ 
cordingly can absorb e/3 dipoles with energetic impunity. 
Incidentally, this is how tunneling fuses the quantum Hall 
fluids together; bulk anyons can then move coherently 
across the trench (leaving behind ‘invisible’ dipoles) even 
though the intervening region supports only electrons. 

To correctly count the degeneracy encoded by the zero 
modes, it is essential to observe that the above quasipar¬ 
ticle operators do not commute, 

[ l p2e/3 5 Vklipole { x )] 7 ^ 0? (13) 

reflecting the fractional statistics of the corresponding 
bulk anyons. One therefore can not simultaneously spec¬ 
ify a unique phase for both (V> 2 e/ 3 ) and (V’dipoie) even 
though the condensates appear in different spatial re¬ 
gions. In a basis where the Pdipoie integers lock to 
specific values the P 2 e /3 integers fluctuate wildly. The 
dipole charges on tunneling-gapped regions correspond¬ 
ingly fluctuate through the allowed low-energy configura¬ 
tions while each paired region realizes degenerate states 
with well-defined total charge of 0, 2e/3, or 4e/3 (mod 
2e). Conversely, specifying p 2e /3 requires Pdipoie to fluc¬ 
tuate, and in this basis the charges on paired regions fluc¬ 
tuate while each tunneling-gapped region takes on one 
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of three degenerate configurations characterized by well- 
defined dipole charges. In short, we can label ground 
states by charge configurations of either the tunneling 
or pairing regions, but not both. Introducing 2N do- 

mains walls then yields a degeneracy of y/3 (neglecting 
Hilbert-space constraints). It follows that the ground- 
state degeneracy per domain wall—also known as the 
quantum dimension—equals y/3. 

The Z 3 parafermion zero-mode operators [denoted ol\^ 
in Fig. [2j[a)] simply cycle the system through the ground- 
state manifold. To gain further intuition for these modes, 
consider adding charge 2e/3 to one end of the trench. 
Doing so within a pairing-gapped region costs energy 
(because of the ensuing dipole moment) and similarly 
for the tunneling regions (due to the net excess charge 
created). Within a domain wall, however, the following 
nontrivial process can uniquely arise: (i) The added 2e/3 
charge propagates chirally until it hits a superconduct¬ 
ing region and Andreev reflects into a — 2e/3 hole on the 
opposite side of the trench. The paired condensate ab¬ 
sorbs the deficit 4e/3 charge, (ii) The —2e/3 hole hops 
back across the trench upon impinging on a tunneling- 
gapped region, whose condensate absorbs the resulting 
change in dipole charge. Figure [2](a) sketches both pro¬ 
cesses. (Hi) Repeating the preceding steps converts the 
hole back into a 2e/3 charge. Parafermion zero-mode 
operators essentially create a uniform superposition of 
these configurations, which costs no energy.^ This phys¬ 
ical picture shows that the zero modes indeed localize to 
domains walls, add fractional charge of 2e/3, and thus in¬ 
herit their nontrivial commutation relations directly from 
the anyons native to the setup—consistent with the in¬ 
tuition laid out at the beginning of this subsection. 


2. Zero-mode detection 

Armed with these insights, one can anticipate several 
sharp experimental signatures of parafermions. Imag¬ 
ine bringing a domain wall to the outer boundary of the 
host v m 2/3 quantum Hall fluid as in Fig. |3][a), thus 
allowing gapless fractionalized edge states to hybridize 
with a parafermion zero mode. In the asymptotic low- 
energy limit, the zero mode mediates ‘perfect Andreev 
conversion’ Incident 2e/3 quasiparticles convert with 
unit probability into — 2e/3 holes at the domain wall, 
with the adjacent superconducting region absorbing the 
deficit 4e/3 charge. Two-terminal transport from the 
edge through the superconductor is thus characterized 
by a quantized zero-bias conductance G = 2 ue 2 /h—twice 
the value obtained when current passes between two nor¬ 
mal leads. Perfect Andreev conversion represents a chiral 
analogue of the well-studied Majorana-induced perfect 
Andreev reflection 53 56 ^ and provides a conceptually sim¬ 
ple fingerprint of parafermion zero modes. This effect 
also offers practical utility, as one can exploit its chiral 
nature to build exotic circuit elements including transis¬ 
tors for fractional charge and flux capacitorsP° I 


A variant of Fig. [3][a) further enables detection of 
the domain wall’s quantum dimension through tunneling 
spectroscopy. In particular, introducing a second domain 
wall near the outer edge couples a pair of zero modes, 
thereby reducing the degeneracy. The zero-bias conduc¬ 
tance peak noted above consequently splits into g finite- 
voltage resonanced^ with g the square of the quantum 
dimension. 

Consider next Fig. !b) where parafermion modes 
oq 5 2 hybridize across a Josephson junction with phases 
4>l/r on the left/right superconductorsF 2 48 49 58 Here 
the parafermions yield an anomalous Josephson current 
mediated not by tunneling of ordinary Cooper pairs, but 
rather fractional 2e/3 charges. This effect is actually 
quite natural given the aforementioned anyon condensa¬ 
tion present in the paired regions. Viewing e l ^ L / R as a 
charge- 2 e Cooper-pair creation operator, we see that the 
anomalous current (in the weak hybridization limit) takes 
the form 

I oc sin ( — (14) 

where the offset phase 0 depends on the precise state of 
the barrier. The Josephson current, in stark contrast to 
ordinary junctions, thus exhibits 67 r periodicity in the 
phase difference </>z, — </>#! This phenomenon general¬ 
izes Kitaev’s 47 r-periodic ‘fractional Josephson effect’ for 
Majorana systems, where supercurrent can arise via sin¬ 
gle electron tunneling ! 5 The present setup, however, of¬ 
fers a potential advantage. Namely, whereas stray un¬ 
paired electrons spoil the unconventional periodicity in 
the Majorana case, the parafermionic Josephson effect 
is immune to environmental electronic noise since finer 
fragments carry the anomalous supercurrent. 

3. Non-Abelian statistics 

One can profitably view the superconducting region of 
the trench in Fig. [2^a) as an extrinsic line defect em¬ 
bedded brute force into the quantum Hall fluid. The 
endpoints—where parafermion zero modes localize— 
represent topological point defects with the key prop¬ 
erty that a bulk anyon encircling one of the ends swaps 
the sign of its charge. [The sign change occurs via an 
Andreev-like process similar to Fig. [3](a) when the anyon 
crosses the paired region.] Such topological defects are 
confined since they clearly incur an energy cost that 
grows linearly with their separation. Nevertheless it is 
valuable to view them as ‘particles’, particularly since 
it is often possible to devise schemes for controllably 
nucleating and mobilizing parafermion zero modes; see, 
e.g., Refs. H81 and 49 The theory of extrinsic defects in 
fractionalized media (whose properties intimately relate 
to symmetries of the topological field theory describing 
the host platform) has been developed in great detail, 
e.g., in Refs. QTJ [42j [48j [49j [5TJ I55H701 Here we will 
simply address the following important question: What 





(a) 


(b) 



FIG. 3. Detection schemes for parafermion zero modes, illustrated using the v — 2/3 setup from Fig. [ 2 ^a). (a) Hybridizing a 
zero mode with the gapless edge states generates ‘perfect Andreev conversion’, (b) Coupling Z 3 parafermions across a Josephson 
junction yields a current (mediated by 2e/3 tunneling) 67 r-periodic in the phase difference 4>l — 4>r- 


kind of ‘particle’ is the topological defect that binds a Z 3 
parafermion zero mode in our v = 2/3 setup? 

Suppose that we introduce a set of 27V well-separated 
topological defects into the quantum Hall fluid, generat- 

ing y/3 ground states (again neglecting Hilbert-space 
constraints). Importantly, this degeneracy does not arise 
from a local property of the defects, consistent with their 
non-integer quantum dimension, but instead reflects the 
degenerate global charge configurations on the supercon¬ 
ducting regions linking parafermion zero modes. Hence 
provided the zero modes maintain their distance no lo¬ 
cal perturbation can split the ground states; they enjoy 
topological protection. Even more interestingly, adiabat- 
ically exchanging pairs of defects noncommutatively ro¬ 
tates the system’s quantum state within the degenerate 
manifold. In other words, the defects binding parafermion 
zero modes effectively realize non-Abelian anyons —a fas¬ 
cinating conclusion given that the host v = 2/3 system 
comprises a prototypical Abelian quantum Hall state. To 
demystify this seemingly paradoxical statement it is use¬ 
ful to think of the v = 2/3 state as descending from a 
parent non-Abelian quantum Hall phase [closely related 
to SU( 2 ) 4 p 5 71 upon a phase transition that confines the 
non-Abelian anyons. By introducing superconducting re¬ 
gions into the fluid, we simply provide the energy nec¬ 
essary to pull apart these confined non-Abelian excita¬ 
tions, whose exotic properties we can still, remarkably, 
harness.^ 

Widespread interest in non-Abelian anyons stems 
partly from their potential utility for inherently fault- 
tolerant topological quantum computing . 1 2 The key idea 
is that the degenerate ground-state manifold non-locally 
encodes topological qubit states that can be controllably 
manipulated through defect braiding by virtue of their 
non-Abelian exchange statistics. Braiding Majorana 
zero modes imparts only 90° qubit rotations, thus en¬ 
abling somewhat restricted fault-tolerant quantum infor¬ 
mation processing.® It is therefore natural to ask whether 
our more exotic parafermion zero modes offer compara¬ 
tive advantages. While in the parafermion case braid¬ 
ing provides additional gates that become trivial with 
Major anas SSI their exchange statistics remains insuffi¬ 
cient for fully universal topological quantum computa¬ 


tion. We note as a teaser that it is possible to leverage 
closely related setups to a chieve bona fide computational 
universality (see Sec. IV), but first we explore promising 
alternative settings for realizing parafermion zero modes. 


B. Quantum Hall bilayers 


Zooming out slightly, we saw in Sec. Ill A 1 that 
parafermion zero modes arise at the interface between 
two ‘incompatibly’ gapped regions of a trench in a frac¬ 
tionalized medium. The Cooper pairing exploited above 
provides one gapping mechanism that differs fundamen¬ 
tally from ordinary electron backscattering, but inter¬ 
estingly is not always necessary. Certain quantum Hall 
phases admit multiple charge-conserving gap-opening 
processes that allow parafermion zero modes to emerge 
even without superconductivity. 

Figure [2jb) depicts an important example consisti ng of 
a quantum Hall bilayer with v m 1/3 filling per layer ESEH 
A trench through the system can acquire a gap via ( 1 ) ‘di¬ 
rect’ tunneling, with electrons backscattering top-to-top 
and bottom-to-bottom and (ii) ‘crossed’ tunneling where 
electrons backscatter by moving between layers. Viewing 
the layer index as a pseudospin, these mechanisms gap 
the charge sector in the same fashion (inevitably, because 
charge is conserved) but produce distinct gapped phases 
in the pseudospin channel. 

Domain walls separating these incompatibly gapped 
phases once again trap Z 3 parafermion zero modes. The 
associated ground-state degeneracy arises because the 
hopping perturbations catalyze different flavors of anyon 
condensation that allow the trench to freely absorb ‘hori¬ 
zontal’ intra-layer and ‘diagonal’ inter-layer e/3 dipoles in 
direct- and crossed-tunneling regions, respectively. Since 
the operators creating these two dipole types do not 
commute—again a reflection of fractional statistics of the 
corresponding bulk anyons—one can label ground states 
by the charge differences (mod 2 e) on either direct- or 
crossed-tunneling regions but not both. The zero-mode 
operators cycle the system through the degenerate mani¬ 
fold by adding ‘vertical’ e/3 dipoles to one side of the 
trench; such dipoles cost finite energy everywhere ex- 











9 


cept inside of the domain walls, which explains why the 
parafermions localize. 

As in the superconducting case discussed previously, 
one can view the endpoints of crossed-tunneling regions 
as confined topological point defects that carry quan¬ 
tum dimension y/3 and exhibit non-Abelian statistics. 
A bulk anyon winding around such a defect ends up 
in the opposite layer from which it began (rather than 
reversing its charge). Given the quantum dimension, 
each additional crossed-tunneling region inserted into the 
^ = l/3 + l/3 bilayer increases the ground-state degener¬ 
acy threefold. Curiously, the degeneracy of a single-layer 
v = 1/3 fluid increases by precisely the same factor when 
system’s genus rises. This is no coincidence—upon rotat¬ 
ing the top layer in Fig. |2|b) by 180° about the trench, 
the crossed-tunneling region locally resembles a puncture 
drilled into ina^ = 1/3 system on a sphere. Barkeshli et 
alP° thus christened the topological defects in this bilayer 
setup ‘genons’. 

Parafermion zero modes bound to ‘genons’ yield trans¬ 
port fingerprints similar to those arising in quantum 
Hall/superconductor heterostructureFor instance, 
hybridizing a zero mode with gapless edge states na¬ 
tive to the bilayer [in a straightforward generalization 
of Fig. [3ja)] generates a perfect-Andreev-reflection ana¬ 
logue where incoming ‘vertical’ dipoles from the edge 
switch polarity with unit probability.^ (The adjacent 
crossed-tunneling condensate readily absorbs the change 
in dipole moment.) More broadly, the similarity between 
the two setups discussed so far is quite remarkable given 
the total absence superconductivity in the present case. 
While the bilayer architecture nevertheless raises nontriv¬ 
ial fabrication challenges given the somewhat unnatural 
crossed-tunneling terms that we invoked, Ref. 731 suggests 
several plausible experimental routes to realizing the cou¬ 
plings (interestingly, not all of which are essential for the 
topological degeneracy). 


C. 2D topological insulator edges 

While we have now seen that fractionalized 2 D systems 
with anyons indeed comprise natural parafermion-zero- 
mode hosts, fractionalization turns out, quite remark¬ 
ably, to be inessential. Two important works TMIM re _ 
cently showed that one can instead trap parafermions us¬ 
ing a much simpler medium: a 2D topological insulator. 
As background, a 2 D topological insulator edge supports 
counterpropagating electron modes described by opera¬ 
tors fn/L that form Kramers partners. These modes nec¬ 
essarily remains gapless provided the system preserves 
both particle-conservation and time-reversal symmetries. 
Conversely, adding superconductivity induces a pairing 
gap while introducing magnetism can gap the edge via 
backscattering. As shown earlier by Fu and Kand^ 
domain walls separating these clearly distinct gapped 
phases bind Majorana zero modes. More interestingly, 
replacing the magnetic backscattering terms with strong 


interactions that spontaneously break time-reversal sym¬ 
metry promotes these Majorana modes to more exotic Z 4 
parafermions! 

We can capture this result using a bosonized frame¬ 
work where the electron operators and density respec¬ 
tively become /r/l ~ e l ^ ±e ) and p = d x 6 /ir. The 
bosonized fields obey [(p(x),0(x f )\ = mQ{x — x' ) with 
Q(x) the Heaviside function and x an edge coordinate. 
Let us examine in this language the gapping mechanisms 
available in the absence of explicit time-reversal break¬ 
ing: (i) Proximity-induced Cooper pairing with strength 
A yields a perturbation 

tfpair = - [ (A/r/l + H.C.) - f Acos(2 ip) (15) 

J X J X 

that gaps the edge by locking ip = 7 rp pa i r with p pa i r €= ^ 
to minimize the energy, (it) Time reversal forbids elas¬ 
tic single-electron backscattering but allows two-particle 
backscattering described by a Hamiltonian 

H hs = - f (Tfyy L f L + H.c.) ~ - f r cos(4$). (16) 

J X J X 

With sufficiently strong density-density repulsion the 
coupling T becomes relevant and gaps the edge by pin¬ 
ning 9 = |pbs where pbs £ ^JHHEDThe system then spon¬ 
taneously breaks time reversal symmetry, as indicated 
by the local magnetization order parameter (/r/l) ~ 
(—l) Pbs . Contrary to the situation with explicit sym¬ 
metry breaking, however, configurations with opposite 
magnetization (i.e., even vs. odd pbs) remain degenerate. 
This feature proves essential in what follows. 

Consider the setup of Fig. [2|c) where regions gapped 
by two-electron backscattering straddle a Cooper-paired 
edge segment proximate to a superconductor l 77 The 
backscattering terms pin 0 = fpj/sT on left/right re¬ 
gions while cp = 7 rp pa i r on the superconducting section. 
Similar to the situation encountered in Sec. |III A H we 
can not uniquely fix all of these integers since 0 and cp do 
not commute. The special physics of this setup becomes 
most apparent in a basis where we specify p 1 ^, leaving 
the p pa i r integers wildly fluctuating. The total charge 
(mod 2 e) on the central paired segment then reads 


Q = e r 

Jxi 


' d x 9 0(x r ) - 0(xi) e t 

dx - = e-= -(Pbs -Pbs), 

7r 7T Z 


(17) 


with xi/ r positions just to the left/right of the super¬ 
conductor. Hence the superconducting part of the edge 
admits four degenerate states labeled by total charges 
0, e/2, e, and 3e/2 (mod 2e). Furthermore, we see that 
raising p£ s or p l hs by one unit—thereby reversing the mag¬ 
netization in an adjacent backscattering region—pumps 
charge ±e/2 into the superconductor. Accordingly, local¬ 
ized Z 4 parafermion zero-mode operators cycle the sys¬ 
tem through the degenerate manifold by adding spin and 
charge e /2 to the domain walls 
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Hybridizing a pair of Z 4 parafermions across a Joseph- 
son junction [similar to Fig. |3jb)] yields an anomalous Sir- 
periodic Josephson current reflecting e/2 transport across 
the barrier. Perhaps the most remarkable feature of this 
setup is that while isolated parafermion zero modes re¬ 
quire strong electron-electron repulsion—a nontrivial sce¬ 
nario to obtain in practice— the Sir-periodic Josephson ef¬ 
fect can survive even with arbitrarily weak interactions 
This beautiful result may be understood by adiabatic 
continuity: Starting from the strongly coupled limit, di¬ 
aling down the barrier’s interaction strength yields only 
smooth crossover behavior due to finite-size effects, hence 
the robustness of the Sir periodicity. 

Since trapping parafermion zero modes does not re¬ 
quire systems with anyons, one might hope that even 
simpler purely ID wires suffice. Several works have ex¬ 
plored this intriguing issue 79 -^^, providing possible routes 
that do, however, require some amount of fine-tuning—in 
accord with the ID classifications discussed at the start 
of Sec. |III| If parafermion zero modes can stably appear 
in strictly ID setups they will likely require symmetry 
protection; finding explicit examples (or proving other¬ 
wise) remains an important open question. 


IV. FIBONACCI ANYONS FROM COUPLED 
PARAFERMIONS 


Braiding topological defects binding parafermion zero 
modes yields denser qubit rotations compared to the Ma- 
jorana case, but not sufficiently dense to enable compu¬ 
tational universality. In search of universal topological 
quantum computing hardware we must therefore look to 
still more exotic non-Abelian excitations. Following the 
spirit of the above analysis, here we discuss a way to 
“engineer” the requisite topological order. 

For inspiration, we recall that potential “intrinsic” 
sources may arise from highly exotic quantum Hall 
phases. The Z 3 Read-Rezayi stated, a natural gener¬ 
alization of the famed Moore-Read statd^, provides a 
notable example. This phase features a gapless edge 
structure consisting of a charge mode together with a 
neutral sector described by a chiral Z 3 parafermion con¬ 
formal field theory (whic h com prises half of the non-chiral 
theory discussed in Sec. IIC). As an intimately related 
consequence the bulk supports non-Abelian ‘Fibonacci 
anyons’ with the following properties: their quantum 
dimension is the golden mean, cp = (1 + \/5)/2; two 
Fibonacci anyons brought together either annihilate or 
beget another Fibonacci anyon; and, most in terest ingly, 
their braiding provides a universal gate set! 86 * 87 . Fi¬ 
bonacci anyons thus constitute a holy grail for topological 
quantum computation. Ultimately the special physics 
of the Read-Rezayi state derives from a peculiar prop¬ 
erty built into the ground state: clustering correlations 
amongst triplets of electrons (as opposed to pairs that 
occur in a superconductor). 

The comparatively simple Abelian quantum Hall sys¬ 


tems exploited to obtain Z3 parafermion zero modes in 
Sec. [HI] can, remarkably, provide close cousins of the 
Read-Rezayi state. As a concrete example, imagine pat¬ 
terning a 2D array of superconducting islands into a 
v = 2/3 state, thereby generating a network of Z 3 zero 
modes that encode a massive ground-state degeneracy 
reminiscent of the situation in a partially filled Lan¬ 
dau level. Judiciously hybridizing adjacent parafermion 
zero modes can then lift the degeneracy in favor of 
a superconducting ‘Fibonacci phase’ hosting a gapless 
chiral Z3 parafermion mode at its boundary and Fi¬ 
bonacci anyons in the interior . 35 This outcome is actu¬ 
ally quite intuitive. As emphasized previously, the su¬ 
perconducting islands impart condensation of charge- 2 e 
Cooper pairs that derive from conglomerates of three 
fractionalized 2e/3 quasiparticles in the quantum Hall 
fluid. Thus the heterostructure effectively leverages pair¬ 
ing to mimic multi-particle clustering—making Read- 
Rezayi-like physics rather natural. The resulting corre¬ 
spondence between the superconducting Fibonacci phase 
and the Read-Rezayi state [summarized in Fig. |4](a)] gen¬ 
eralizes Read and Green’s^ classic work demonstrating 
the common non-Abelian content of a p + ip supercon¬ 
ductor and the Moore-Read state. 

Analytical and numerical analyses substantiate the 
preceding physical argument. One powerful method as¬ 
sembles the 2 D Fibonacci phase from coupled chains sim¬ 
ilar to Teo and Kane’s seminal workP^that obtained non- 
Abelian quantum Hall phases from Luttinger liquid ar¬ 
rays. Here, however, we utilize the parafermion chains 
described in Sec. [TT| instead of Luttinger liquids. Namely, 
we consider Z3 parafermion operators arranged on a tri¬ 
angular lattice with nearest-neighbor couplings £ 1 ,£ 2 , and 
£3 sketched in Fig.^b)^. The setup could arise either in 
a v = 2/3 state with superconductivity or a z/ = 1/3+1/3 
bilayer with crossed-tunneling defects, among others. 

In the extreme limit t\^ = 0 the problem reduces to 
decoupled Potts chains of Sec. |II A[ with J = / = £3 and 
cf) = 0 = 0. As Fig. [ljc) illustrates each chain therefore 
realizes a non-chiral Z3 parafermion critical point. The 
system’s gaplessness here is a great virtue since even ex¬ 
ceedingly weak interchain couplings t\^ can profoundly 
impact the physics. Crucially, the relation between lat¬ 
tice parafermion operators and fields in the correspond¬ 
ing conformal field theory is knowrP^ recall Eqs. ^ 
and (10). This relation allows a systematic analysis of 
the weakly-coupled-chain limit that ‘filters out’ irrele¬ 
vant lattice-scale details by projecting onto low-energy 
fields describing long-wavelength physics. (Reference l89l 
pursues a similar strategy, though we caution that their 
construction utilizes unphysical interchain couplings; see 
Ref. |24jfor further discussion.) 

An especially desirable limit arises when t\ = £ 2 = £, 
where the most relevant interchain perturbation becomes 
simpl^l35|83J 


5H ~ 





(18) 
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FIG. 4. (a) Correspondence between the Z 3 Read-Rezayi state and the ‘Fibonacci phase’ formed in a is = 2/3 quantum 

Hall/superconductor heterostructure. (b) Triangular-lattice parafermion model (left) that realizes the Fibonacci phase in the 
weakly-coupled-chain limit (right), (c) Phase diagram of the setup in (b) obtained from the DMRG analysis of Ref. l83l 


with ^ L/R, y parafermion fields for chain y. The wavy 
lines on the right side of Fig. Qb) designate the pattern 
of couplings, which completely gap the interior but leave 
behind ‘unpaired’ chiral Z 3 parafermion modes at the 
boundary. Such a perturbation of the parafermion criti¬ 
cal point in purely ID field theory (i.e., if we instead cou¬ 
pled 'i/jr to VfL in the same chain), results in a Fibonacci- 
type spectrum, as known exactly from integrability 90 i In 
our coupled-chain system, the perturbation correspondly 
triggers precisely the Fibonacci phase anticipated above. 
While this analytic approach breaks down at stronger in¬ 
terchain couplings, density matrix renormalization group 
(DMRG) simulations allow a controlled exploration of 
the broader parameter spaced Figure [4|c) sketches the 
numerically determined phase diagram; fortunately the 
Fibonacci phase enjoys a wide stability window extend¬ 
ing from weakly coupled chains to the isotropic triangular 
lattice limit (t\ = £2 = £ 3 ) and beyond! 

It is quite encouraging that a non-Abelian state suit¬ 
able for universal topological quantum computation can 
in principle emerge upon combining well understood com¬ 
ponents such as a v = 2/3 state and conventional su¬ 
perconductors. Similar phases have also been predicted 
in various related settings including uniform quantum 
Hall bilayers 9 ^, quantum wire 92 and spin chain 93 ar¬ 
rays, and even local bosonic model^. We also note 
that an interesting earlier study, Ref. 95. explored other 
coupled 2D parafermion lattice models and captured an 
Abelian topologically ordered phase that generalizes the 
toric code. These developments raise a number of in¬ 
triguing questions: How can the architectures reviewed 
here be further distilled to ease practical implementa¬ 


tion? What probes enable unambiguous detection of Fi¬ 
bonacci anyons in such setups? How might they be ma¬ 
nipulated to reveal their coveted braiding statistics? And 
can one pursue similar philosophies to obtain perhaps 
even more exotic non-Abelian phases? 


V. OUTLOOK 

Parafermions boast a more than three-decade history 
in condensed matter, with numerous now-classic applica¬ 
tions ranging from criticality in Potts chains to the non- 
Abelian Read-Rezayi quantum Hall phases. Despite this 
lengthy history, understanding of parafermions—both on 
a formal and practical level—has advanced significantly 
during the past few years. Much of this progress can 
be viewed as extending the wildly influential Ising-model 
physics (in its myriad of incarnations) into more exotic 
territory. 

On the formal end we saw in Sec. [II] that para¬ 
fermion chains exhibit a topological phase supporting 
parafermion zero modes that comprise natural Majorana 
generalizations. While appealing parallels with the sim¬ 
pler Ising/Majorana chain do exist, we encountered sub¬ 
tle surprises owing to the inherently strongly interacting 
nature of parafermions. As one intriguing example, even 
deep within the topological phase strong parafermion 
zero modes (which guarantee threefold degeneracy of all 
eigenstates) can abruptly transform into weak zero modes 
due to an unusual form of quantum criticality operating 
solely within the excited states. The effect of disorder on 
this transition raises an interesting issue that warrants 
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further exploration; in particular, following Refs. [3T1 and 
33 it is tempting to anticipate a many-body localized 
phase that resurrects strong zero modes even when they 
are absent in the clean system. 

Another recent formal advance established the pre¬ 
cise correspondence between lattice parafermion opera¬ 
tors and continuum parafermion fields emerging in the 
critical chain. These relations not only clarify certain 
aspects of the conformal fields theory, but also underlie 
controlled analytical techniques that attack 2D systems 
built from parafermion chains. Yet another example is 
the identification of a ‘free-parafermion’ limit obtained in 
a deformed non-Hermitian parafermion chain model . 30 
Might it be possible to bootstrap off of this limit to 
controllably attack strongly interacting parafermion sys¬ 
tems? Or can one perhaps design implementations of the 
free-parafermion model itself? 

By now many proposed experimental realizations of 
parafermion zero modes exist. We focused in Sec. |III| on 
three classes of systems that appear especially promis¬ 
ing in that they rely entirely on well-understood phases 
of matter: Abelian quantum Hall/superconductor hy¬ 
brids, simple quantum Hall bilayers with ‘genon’ defects, 
and strongly interacting 2D topological insulator edges 
coupled to superconductivity. The first class perhaps 
appears a somewhat awkward marriage given the high 
fields usually required to observe fractional quantum Hall 
physics. We stress, however, that very recent experi¬ 
ments successfully induced a proximity effect in GaAs 
quantum wells using high-critical-field (> 10T) NbN su¬ 
perconducting contacts.^ This development bodes well 
for the prospects of exploring the interplay between quan¬ 
tum Hall physics and superconductivity—a largely un¬ 
tapped experimental arena ripe for discovery. Bilayer 
‘genon’ proposals require controlled manipulation of edge 
states to selectively hybridize different pairs of coun- 
terpropagating modes, which although nontrivial seems 
achievable with presently available gating technology.^ 
Such systems carry the added bonus that detecting 
parafermion zero modes can also reveal (albeit indirectly) 
the topological ground state degeneracy of quantum Hall 
phases on nontrivial-genus manifolds . 57 Interacting 2D 
topological insulator edges may offer the most immedi¬ 
ate prospects of detecting signatures of parafermion zero 
modes. Two virtues are particularly noteworthy here: 
First, experiments have al ready achieved the required 
coupling to superconductors pZIMI and second, the anoma¬ 
lous 87 r-periodic Josephson effect that parafermions me¬ 
diate survives even in the weakly interacting limit ! 74 


One can view arrays of coupled parafermion zero 
modes in such setups as new low-energy degrees of free¬ 
dom from which even more exotic higher-dimensional 
topological phases can emerge. Our understanding of 
the possibilities is certainly far from complete, though 
initial studies have already revealed exceptionally rich 
physics. Notably, networks of parafermion modes arising 
in the above quantum-Hall-based platforms can yield a 
‘Fibonacci phase’ that represents an analogue of the cov¬ 
eted Z 3 Read-Rezayi state with all but the most inter¬ 
esting non-Abelian Fibonacci anyons—which carry uni¬ 
versal braid statistics—distilled away. This insight pro¬ 
vides proof of principle that well understood components 
can combine to yield universal topological quantum com¬ 
puting hardware. Simplifying the requisite architectures 
poses an important outstanding challenge. One promis¬ 
ing avenue is to wield numerics to examine the phase 
diagrams of uniform quantum Hall/superconductor hy¬ 
brids and quantum Hall bilayers in search of related non- 
Abelian phases; studies in the latter setting have in fact 
recently been carried out" 101 . 

While experimental prospects and tantalizing long¬ 
term technological applications are certainly exciting, in 
our view it is worth emphasizing the inherent theoretical 
beauty underlying the advances in parafermions surveyed 
here. The physics that has so far emerged intimately 
links seemingly disparate areas across condensed-matter 
physics to an extent that rivals the reach of the Ising 
model. As this story continues to evolve, undoubtedly 
many more fundamental insights will be gained in this 
fascinating field. 
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